Rotating Optical Soliton Clusters 
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We introduce the concept of soliton clusters - multi-soliton bound states in a homogeneous bulk op- 
tical medium, and reveal a key physical mechanism for their stabilization associated with a staircase- 
like phase distribution that induces a net angular momentum and leads to cluster rotation. The 
ringlike soliton clusters provide a nontrivial generalization of the concepts of two-soliton spiraling, 
optical vortex solitons, and necklace- type optical beams. 



PACS numbers: 42.65.Tg 

Recent progress in generating spatial optical solitons in 
nonlinear bulk media opens the possibility to study truly 
two-dimensional self-trapping of light and interaction of 
multi-dimensional solitary waves Q. Robust nature of 
spatial solitons that they display in interactions allows 
us to draw a formal analogy with atomic physics and 
treat spatial solitons as "atoms of light". Our motiva- 
tion here is to find out whether more complex objects, 
viewed as "atom clusters" , can be constructed from a cer- 
tain number of simple solitons - "atoms" . In this Letter, 
we describe, for the first time to our knowledge, the ba- 
sic principles for constructing the so-called "soliton clus- 
ters" , ringlike multi-soliton bound states in a bulk media. 

First, in order to discuss the formation of multi-soliton 
bound states in a homogeneous bulk medium, we should 
recall the physics of the coherent interaction of two spa- 
tial solitons. It is well known |jj that such an interaction 
depends crucially on the relative soliton phase, say 0, so 
that two solitons attract each other for 6 = 0, and they 
repel each other for = it. For the intermediate values 
of the soliton phase, < < 7r, the solitons undergo an 
energy exchange and display inelastic interaction. As a 
result, no stationary bound states of two coherently in- 
teracting solitons are possible in a bulk medium. 

Soliton spiraling was suggested theoretically |2[ and 
observed experimentally ||] as a possible scenario for 
a dynamical two-soliton bound state formed when two 
solitons are launched with initially twisted trajectories. 
However, it was clarified later || that the experimental 
observation of the soliton spiraling was possible due to 
an effectively vectorial beam interaction. As a matter of 
fact, the soliton spiraling reported in Refs. (|^] is asso- 
ciated with large-amplitude oscillations of a dipole-mode 
vector state ||] generated by the interaction of two ini- 
tially mutually incoherent optical beams. 

In spite of the fact that no bound states exist for two 
coherently interacting scalar solitons in a bulk medium, 
in this Letter we demonstrate that such bound states (or 
ringlike clusters) are indeed possible for larger number of 
solitons, namely for N > 4. The main reason for the ex- 
istence of such multi-soliton states can be explained with 
the help of simple physics. Indeed, let us analyze possible 



stationary configurations of TV coherently interacting soli- 
tons in (2+l)-dimensions. The only finite-energy struc- 
tures that would balance out the phase-sensitive coherent 
interaction of the neighboring solitons should possess a 
ringlike geometry. However, a ringlike configuration of 
N solitons will be radially unstable due to an effective 
tension induced by bending of the soliton array. Thus, a 
ring of TV solitons will collapse, if the mutual interaction 
between the neighboring solitons is attractive, or other- 
wise expand, resembling the expansion of the necklace 
beams ||. Nevertheless, a simple physical mechanism 
will provide stabilization of the ringlike configuration of 
N solitons, if we introduce an additional phase on the 
scalar field that twists by 2i\m along the soliton ring. 
This phase introduces an effective centrifugal force that 
can balance out the tension effect and stabilize the ring- 
like soliton cluster. Due to a net angular momentum 
induced by such a phase distribution, such soliton clus- 
ters rotate with an angular velocity which depends on 
the number of solitons and phase charge m. 

To build up the theory of the soliton clusters, we con- 
sider a coherent superposition of N solitons with the en- 
velopes G n (x,y,z), n = 1,2. .TV, propagating in a self- 
focusing homogeneous bulk medium. The equation for 
the slowly varying field envelope E = ^2 Gn can be writ- 
ten in the form of the nonlinear Schrodinger equation, 
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where is the transverse Laplacian and z is the prop- 
agation distance measured in the units of the diffraction 
length. Function f(I) describes the nonlinear properties 
of an optical medium, and it is assumed to depend on 
the total beam intensity, / = \E\ 2 . 

General features of the dynamical system (0) are de- 
termined by its conservation laws, or integrals of mo- 
tion: the beam power, P = J \ E\ 2 dr, the linear momen- 
tum, L = Im j E*\7Edr, and the angular momentum, 
Me z — Im J £*(rx VE)dr. For a ring of identical weakly 
overlapping solitons launched in parallel, we can calculate 
the integrals of motion employing a Gaussian ansatz for 
a single beam G n , 
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G n = A exp 
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where r n = {x n ;y n } defines the position of the soliton 
center, and a n is the phase of the n-th beam. Then, the 
integrals of motion take the form: 
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P = ira 2 A 2 e Ynk cosO nkj 

n,k=l 
N 

L = |-4 2 e Ynk (r n -r k )smO nk , 

n,k=l 
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M = ttA 2 e Ynk \r n xr k \s'me nk , 
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where Y nk = -|r n - r k \ 2 /4a 2 and nk = a n - a k . We as- 
sume that the beams are arranged in a ring-shaped array: 
r n = {Rcoscp n ; Rsincpn} with cp n = 2im/N, for which 
we find Y nk = -(R/a) 2 sm 2 [7r(n - k)/N). 

First of all, analyzing many-soliton clusters, we remove 
the centre-of-the mass motion and take L = 0. Applying 
this constraint to Eqs. (||), we find the conditions for the 
soliton phases, Qfi+ n — c\i = a k + n — a kl which are easy 
to satisfy provided the phase a n has a linear dependence 
on n, i.e. a n = 9n, where is the relative phase between 
two neighboring solitons in the ring. Then, we employ 
the phase periodicity condition taken in the form a n+ n = 
a n + 27rm, and find: 

(4) 

In terms of the classical fields, Eq. (|J) gives the con- 
dition of the vanishing energy flow L = 0, because the 
linear momentum L = J jdr can be presented through 
the local current j = kn(E*VE). Therefore, Eq. (@) de- 
termines a nontrivial phase distribution for the effectively 
elastic soliton interaction in the ring. In particular, for 
the well-known case of two solitons (N = 2), this condi- 
tion gives only two states with the zero energy exchange, 
when m is even (0 = 0, mutual attraction) and when m 
is odd (6 = 7r, mutual repulsion) |l|]. 

For a given TV > 2, the condition ([|) predicts the 
existence of a discrete set of allowed stationary states 
corresponding to the values = 0^ with m = 
0, ±1, . . . , ±(N - 1). Here two states 0(±l m l) differ only 
by the sign (direction) of the angular momentum, sim- 
ilar to the case of vortex solitons. Moreover, for any 
positive (negative) m + within the domain tt < \6\ < 2tt, 
one can find the corresponding negative (positive) value 
ra_ within the domain < \0\ < 7r, so that both m + 
and vrt- describe the same cluster. For example, in the 
case N = 3, three states with zero energy exchange are 
possible: 0<°) = 0, 6>W = 2tt/3, and 0^ = 4tt/3, and the 
correspondence is O^ 1 ^ 0^ 2 \ Therefore, it is useful 
to introduce the main value of in the domain < < 7r, 
keeping in mind that all allowed states inside the domain 



< < 7r are degenerated with respect to the sign of 
the angular momentum. The absolute value of the an- 
gular momentum vanishes at both ends of this domain, 
when m = 0, for any A/", and when m = N/2, for even N. 
The number m determines the full phase twist around 
the ring, and it plays a role of the topological charge of 
a phase dislocation associated with the ring. 

In order to demonstrate the basic properties of the soli- 
ton clusters for a particular example, we select the well- 
known saturable nonlinear Kerr medium with F(I) = 
1(1 + s/) _1 , where s is a saturation parameter. This 
model supports stable (2+l)-dimensional solitons. First, 
we apply the variational technique to find the parame- 
ters of a single soliton described by the ansatz (||), and 
find A = 3.604 and a — 1.623 for s = 0.5. Then, substi- 
tuting Eq. (||) into the system Hamiltonian, we calculate 
the effective interaction potential U (R) = H(R)/\H(oo)\, 
where H is the system Hamiltonian, 
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As a result, for any N we find three distinct types of 
the interaction potential U(R), shown in Fig. 1 for the 
particular case N = 5. Only one of them has a local min- 
imum at finite R which indicates the cluster stabilization 
against collapse or expansion. 
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FIG. 1. Examples of the effective potential U(R) for a cir- 
cular array of N = 5 solitons. Corresponding values of the 
topological charge m are shown near the curves. Dynamically 
stable bound state is possible for m = 1 only. 

To verify the predictions of our effective-particle ap- 
proach, we perform a series of simulations of different N- 
soliton rings, using the fast-Fourier-transform split-step 
numerical algorithm and monitoring the conservation of 
the integrals of motion. Alongside with non stationary 
behavior, such as breathing and radiation emitting, we 
find the clusters dynamics in excellent agreement with our 
analysis. According to Fig. 1, the effective potential is 
always attractive for m = 0, and thus the ring of N = 5 
in-phase solitons should exhibit oscillations and, possi- 
ble, soliton fusion. Indeed, such a dynamics is observed 
in Fig. 2(a). Although the oscillations of the ring are 
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well described by the effective potential U(R), the ring 
dynamics is more complicated. Another scenario of the 
mutual soliton interaction corresponds to the repulsive 
potential [shown, e.g., in Fig. 1 for the case m = 2 and 
6 = 47r/5]. In the numerical simulations corresponding 
to this case, the ringlike soliton array expands with the 
slowing down rotation, as is shown in Fig. 2(d). 
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FIG. 2. Different regimes of the interaction of N — 5 soli- 
tons: (a) m = 0, collapse and fusion through oscillations; (b) 
m = 1, a stationary bound state with Rq = 3; (c) m = 1, an 
excited bound state with the oscillation period z pe riod = 22; 
(d) m = 2, the soliton repulsion. The initial radius of the ring 
in the cases (a), (c), and (d) is R = 5. 

Evolution of the stationary bound state that corre- 
sponds to a minimum of the effective potential U (R) in 
Fig. 1 (for m = 1) is shown in Fig. 2(b). Here the an- 
gular momentum is nonzero, and it produces a repulsive 
centrifugal force that balances out the soliton attraction. 
The effective potential predicts the stationary state at 
Ro = 3 with a good accuracy, and the cluster does not 
change its form while rotating during the propagation. 
To continue the analogy between the soliton cluster and 
a rigid body, we calculate the cluster's moment of inertia, 
X, and its angular velocity, fl: 



E\ 2 r 2 dr, tt = M/T. 



(5) 



For the case shown in Fig. 2(b), the numerically obtained 
value of the angular velocity is £l nU m — tt/20 = 0.157, 
while the formula (||) gives the value Q = 0.154. We 
also perform the numerical simulations of the "excited" 
clusters, as is shown in Fig. 2(c), and observed oscilla- 
tions near the equilibrium state. Such a vibrational state 



of the "TV-soliton molecule" demonstrates the dynamical 
radial stability of the bound state in agreement with the 
effective-particle approximation. 

Our analysis is valid for any TV, and it allows us to 
classify all possible scenarios of the soliton interaction 
in terms of the phase jump between the neighboring 
solitons in the array. Indeed, for = 0, the ring of TV 
solitons collapses through several oscillations. If the main 
value of belongs to the segment < 6 < 7r/2, the 
interaction between solitons is attractive, the value of the 
induced angular momentum is finite, and there exists a 
rotating bound state of TV solitons. However, if belongs 
to the segment tt/2 < 6 < 7r, the soliton interaction is 
repulsive and the soliton ring expands with or without 
(0 = 7r) rotation, similar to the necklace- type beams ||. 
For example, in the case TV = 3 two stationary states are 
possible, = and = 27r/3, and there exist no bound 
states. For TV = 4 and m = 1, the value of is tt/2 and 
a cluster is indeed possible, as is shown in Fig. 3. 
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FIG. 3. Rotating cluster of N = 4 solitons. The parameters 
are Ro = 3.8 and Q = 0.042. Phase images are scaled from 
— 7r (black) to +7T (white). 

Together with the intensity of the four-soliton cluster, 
in Fig. 3 we show the phase distribution for the distances 
up to 60 diffraction lengths. The initial staircase-like 
phase in the ring preserves its shape, and it is a nonlinear 
function of the polar angle (p, similar to the phase of the 
necklace-ring vector solitons with a fractional spin [ff|. 
Note, that the phase of such a state can be described as 
a zeroth- order term in the expansion of the vortex phase 
near the n-th soliton center: 
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Calculating the minimum of the potential U (R) that cor- 
responds to the soliton cluster, we find that for given m 
and N ^> 1, the stationary cluster approaches a vortex 
soliton of the charge m. For example, the equilibrium ra- 
dius Rq for the clusters with m = 1 is Ro = 3.8 for N = 4, 
and for N > 4 it approaches the corresponding vortex ra- 
dius Ro = 3. Furthermore, for m = 2, the soliton bound 
states are possible only if = 4tt/N < tt/2. This gives 
the condition N > 8, and for N > 9 we find Ro = 5 
which is close to the radius of a double-charged vortex. 
Thus, the ringlike soliton cluster can be considered as a 
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nontrivial "discrete" generalization of the optical vortex 
soliton [||. Clusters are generally metastable, and they 
experience the symmetry-breaking instability. However, 
they can propagate for many tenths of the diffraction 
length, being also asymptotically stable in some types of 
nonlinear media similar to the vortex solitons 

In Figs. 4(a,b,d) we show some examples of the excited 
rotating clusters with different number of solitons N and 
initial radius R = N, while in Fig. 4(c) we present the 
stationary rotating ringlike cluster of N = 8 solitons. 
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FIG. 5. Examples of the vector soliton clusters. 

In conclusion, we have revealed a key physical mecha- 
nism for stabilizing multi-soliton bound states in a bulk 
medium in the form of rotating ringlike clusters. Such 
soliton clusters can be considered as a nontrivial general- 
ization of the important concepts of the two-soliton spi- 
raling, optical vortex solitons, and necklace scalar beams, 
and they provide an example of the next generation of 
multiple soliton-based systems operating entirely with 
light. We believe that the basic ideas presented in this 
Letter will be useful for other applications, such as the 
beam dynamics in plasmas and the Skyrme model of 
a classical field theory [Q, and they can be also gener- 
alized to the inhomogeneous systems (such as the Bose- 
Einstein condensates in a trap ]l4|]). 

The authors are indebted to M. Segev and M. Soljacic 
for encouraging comments and suggestions, and E.A. Os- 
trovskaya for a help with numerical simulations and crit- 
ical reading of this manuscript. 



FIG. 4. Examples of the rotating soliton clusters: (a) 
N = 6 and m = 1; (b) N = 7 and m = 1; (c) N = 8 and 
m = 2, here we show the exact bound state with Ro = 5.5 
and the angular velocity Q = 0.092; (d) N = 9 and m = 2. 

We should stress that the staircase-like phase distribu- 
tion is a distinctive feature of the soliton cluster. The 
evolution of the soliton ring with the linear phases, i.e. 
those presented by the first-order terms of Eq. (||), can 
be associated with complex deformation of the vortex. 

More interesting structures are found for the vector 
fields and incoherent interaction of solitons. For exam- 
ple, three coherently interacting solitons cannot form a 
bound state. However, adding a single "atom" to the 
incoherently coupled additional component E\ [see Fig. 
5(a)] leads to mutual trapping of all beams. Here, the 
incoherent attraction balances out both the coherent re- 
pulsion, as in the case of multipole vector solitons [ [lo| , 
and the centrifugal force induced by a net angular mo- 
mentum. This leads to the cluster rotation similar to 
the two-lobe rotating "propeller" soliton ]ll|. In Fig. 
5(b), we show the five- lobe analogue of the necklace-ring 
vector solitons recently discussed in Ref. Q. Unlike the 
necklace-ring solitons, the vector cluster rotates and un- 
dergoes internal oscillations as it propagates. 
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